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Abstract. We consider the diffusion scaling limit of the vicious walker model that is a system of nonintersecting 
random walks. We prove a functional central limit theorem for the model and derive two types of nonintersecting 
Brownian motions, in which the nonintersecting condition is imposed in a finite time interval (0, T] for the first 
type and in an infinite time interval (0, oo) for the second type, respectively. The limit process of the first type is 
a temporally inhomogeneous diffusion, and that of the second type is a temporally homogeneous diffusion that is 
identified with a Dyson's model of Brownian motions studied in the random matrix theory. We show that these two 
types of processes are related to each other by a multi-dimensional generalization of Imhof 's relation, whose original 
form relates the Brownian meander and the three-dimensional Bessel process. We also study the vicious walkers with 
wall restriction and prove a functional central limit theorem in the diffusion scaling limit. 
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1 Introduction 

The system of one-dimensional symmetric simple random walks, in which none of walkers have met others in 
a given time period, is called the vicious walker model. ( See Fisher's paper [8], in which it was introduced 
as a model of statistical mechanics.) The purpose of this paper is to study the scaling limit of vicious walkers 
as a stochastic process. Since each random walk tends to a Brownian motion in the diffusion scaling limit, 
an interacting system of N Brownian motions will be constructed as the scaling limit of vicious walkers with 
an arbitrary finite number of walkers N. We show that a functional central limit theorem for vicious walkers 
holds and the limit process X(t) = {Xi{t), X2{t), . . . ,Xf^{t)) is a temporally inhomogeneous diffusion, that 
is, its transition probability depends on the time interval (0, T] in which the nonintersecting condition is 
imposed. We claim that when N = 2, the process {X2{t) — Xi{t))/\/2 is a one-dimensional Brownian 
motion conditioned to stay positive during a finite time interval (0, T], which is called a Brownian meander 
in [35, 40]. 

We also study the temporally homogeneous diffusion process Y(t) = {Yi{t),Y2{t), . . . ,YN{t)), which is 
obtained from the previous process by taking T — s- oo. The process Y(t) is the Doob ti-transform [5] of the 
absorbing Brownian motion in a Weyl chamber with harmonic function /ijv(x) — Y[i<i<j<Ni-'^J ~ ^"^'^ 
can be regarded as a system of Brownian motions with the drift terms acting as the repulsive two-body 
forces proportional to the inverse of distances between particles [11]. In other words, if we set T — > oo, the 
scaling limit of vicious walkers can realize a Dyson's Brownian motion model studied in the random matrix 
theory [7, 26]. We show the following relation between the processes X(i) and Y(t): 

F(X(-) e dw) ^ cn{T)P{Y{-) e dw) ^ 



h{w{T))' 

where cn(T) is the normalization constant. This equality is a generalization of Imhof 's formula, which 
relates the Brownian meander and the three-dimensional Bessel process [13]. 

The Gaussian ensembles of random matrices can be regarded as the thermodynamical equilibrium of 
Coulomb gas system and that is the reason why Dyson introduced a one-dimensional model of interacting 
Brownian particles with (two-dimensional) Coulomb repulsive potentials [7, 26]. Similar relations between 
our processes and random matrix ensembles can be seen. The distribution of Y(t) is described by using 
the probability density of eigenvalues of random matrices in the Gaussian unitary ensemble (GUE) with 
variance t. Pandey and Mehta [27, 33] introduced a Gaussian ensemble of Hermitian matrices depending on 
a parameter a S [0, 1]. When a = 0, the ensemble is the Gaussian orthogonal ensemble (GOE), and when 
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a = 1, it is the GUE. We will find that the probability density function of ^jT/{t{2T — t)}X(t) coincides 
with that of eigenvalues of matricc;s in the Pandey-Mehta ensemble with q: = ^(T - t)/T. 

We also study vicious walkers with wall restriction and prove the functional central limit theorem in the 
diffusion scaling limit. In this case the obtained temporally inhomogeneous diffusion process is a system of 
nonintersecting Brownian meanders, which is related to the nonstandard classes of random matrices [20, 21]. 



2 Statement of Results 

2.1 Vicious walkers without wall restriction 

Let ({Sj}j>o, -P^) be the A''-dimensional Markov chain starting from z = (2:1, Z2, . . . , zat), such that the 
coordinates S'^ , k = 1,2, ... ,N, are independent simple random walks on Z. We always take the starting 
point z from the set 



{z = (zi, Z2, . . . , zn) G (2Z)^ ; zk+i - Zfc G 2Z+, k ^ 1, . . . , N - 1}, 



(2.1) 



where Z+ is the set of positive integers. Now we consider the condition that any of walkers does not meet 
other walkers up to time m, i.e. 

Sj <S] <■■■< Sf, 0<j<m. (2.2) 

We dente by the conditional probabiHty of under the event = {Sj < Sj < ■ ■ ■ < < j < m}. 
The process ({Sj}j>o, Q^) is called the vicious walkers (up to time to) (see Fisher [8]). 

m=10 
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Figure 1: An example of vicious walks without wall restriction. 

For T > and z e Z^, we consider probability measures n^j', L > 1, on the space of continuous paths 
C([0,T]^R^) defined by 



(2.3) 



where S{t), t > 0, is the interpolation of the random walk S^-, j = 0, 1, 2, . . .. We study the limit distribution 
of the probability /i^ y, i — > 00. 

We put R;^ = {x G : xi < X2 < ■ ■ ■ < xpf}, which is called the Weyl chamber [9]. By virtue of 
the Karlin-McGregor formula [15, 16], the transition density function /Af(t,y|x) of the absorbing Brownian 
motion in R^ and the probability j\/jv(i, x) that the Brownian motion started at x G R< does not hit the 
boundary of R< up to time t > are given by 



/jv(t,y|x) 



det 

l<i,j<N 



((27ri)-'/' e-(-^-''')'/2*) , x,y G R^, 



(2.4) 
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and 



^^N{t,yi)= [ rfy/jv(t,y|x). (2.5) 

For an even integer n and an antisymmetric n x n matrix A = {aij) we put 

Ph<i<j<n{aij) = (•)^^/2)! y^^ff^('^)"<^(l)o-(2)^o-(3)<^(4) • • • ('■a(n-l)a(n), (2-6) 

where the summation is extended over all permutations cr of (1, 2, . . . , n) with restriction a{2k — 1) < <j{2k), 
A; = 1, 2, ... , n/2. This expression is known as the Pfaffian (see Stembridge [39]). Then we have the following 
lemma, which is a consequence of the identity given by de Bruijn [4] as shown in Section 4. 

Lemma 2.1 For t > 0, x e R^, 

\ Ph<i<j<N+iFij{t,x), ifN=odd, ^ ' 

where 



1, ifl<i<N,j = N+l, (2.8) 

-1, ifi = N + l,l<j<N, 

0, ifi = N+l,j = N + l, 



and *(u) = (2/V7r) j^e'^^dv. 
We put 

/iiv(x)= n {xj-Xi). (2.9) 

l<i<j<JV 

The first main result is the following theorem. 

Theorem 2.2 (i) For any fixed z e and T > 0, as L ^ oo, Mz,,t(') converges weakly to the law of 
the temporally inhomogeneous diffusion process X(f) = {X\{t),X2{t), . . . ,Xj^{t)),t G [0, T], with transition 
density g'^{s,x,t,y); 

5^(0,0,t,y) = cjvT^(^-^)/V^'/^exp|-M!|ft^(y)A/V(T-t,y), (2.10) 

.^(.,x,,y)= ^-^^-;y-(^-^'^) , (2.11) 

/or < s < f < T, x,y e R^, w/iere cjv = 2-^/Vnf=i r(j/2). 
(ii) T/ie diffusion process X(f) solves the following equation: 

Xi{t) = Bi{t)+ [\fis,X{s))ds, t€[0,T], i = l,2,...,N, (2.12) 
Jo 

where Bi{t), i = 1,2, . . . , N, are independent one- dimensional Brownian motions and 

d 

bf{t, x) = ^ lnMv(T - t, x), i=l,2,...,N. 
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Although y(-) is the probability measure defined on C([0, T] — > R^), it can be regarded as that on 
C{[0, oo) — > R^) concentrated on the set 

{w e C([0, oo) ^ R^) : w{t) = w{T), t > T}. 

Next we consider the case that T = Tl goes to infinity as L — > oo. 

Corollary 2.3 (i) Let Tl he an increasing function of L with ^ oo as L ^ oo. For any fixed z G Z;^, 
as L ^ 00, ij,ltl(') converges weakly to the law of the temporally homogeneous diffusion process Y{t) = 
{Yi{t),Y2{t), . . . ,YN{t)),t G [0,00), with transition density pN{s,x,t,y); 

Pn{0, 0, t, y) = 4^-^'/^ ^""P {-^} ''^(y^'' (^-l^) 

PN{s,x,t,y) = \ /jv(t- s,y|x)/ijv(y), (2.14) 

for < s <t <oo, x,y G R^, luhere = (27r)-^/V T{j). 

(ii) The diffusion process Y{t) solves the equation of Dyson's Brownian motion model [7] : 

Yi{t)=Bi{t)+ I Yi )-yf ^ = 1,2,...,A^. (2.15) 

2.2 Vicious walkers with wall restriction 

In this subsection, we impose the condition 

Sj > 0, l<j<m, (2.16) 
in addition to (2.2) and take the starting point z from the set 

= {z = (zi, Z2, ...,ZN)e {2Nf ; Zk+i - G 2Z+, fc = 1, . . . , iV - 1}, (2.17) 

where N is the set of non-negative integers. That is, there assumed to be a wall at the origin and all walkers 
can walk only in the region [0,oo). We dente by the conditional probability of under the event 
Am — {0 < Sj < S'j < ■ ■ ■ < Sj^ , < j < m}. The process ({Sj}j>o, Q^) is regarded as the vicious walkers 
with wall restriction (up to time m) [22]. 

For T > and z e N^, we consider probability measures -p, i > 1, on C([0,T] — > R") defined by 

m1,t(-) = Qhr {^^(Lh) e . (2.18) 

We study the limit distribution of the probability juf^ j., L — > 00. 

We put R+< = {x e R^ : < a;i < a;2 < • • • < x^}. Then the transition density function /jv(f, y|x) of 
the absorbing Brownian motion in R+;^ and the probability A/jv(i, x) that the Brownian motion started at 
X e R+< does not hit the boundary of R+< up to time t>0 are given by 

Mt, y|x) = det ({2nt)-'/^ (e-(x.-!/.)V2t _ ^-{^,+mr/2tA , x, y e R+^, (2.19) 

l<i,j<N \ / 

and 

A?'jv(t,x)= / dy/Ar(t,y|x). (2.20) 

Then we have the following lemma, which is proved in Section 4 as a consequence of the identity given by 
de Bruijn [4]. 
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m=10 
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Figure 2: An example of vicious walks with wall restriction. 



Lemma 2.4 For t> 0, x e R+< , 



where 



and 



We put 



J^N{t,x) = 



Ph<t<]<NFij{t,x.), if N = even, 
Pii<i<j<N+iFij{t,x.), if N = odd, 



$(U1,U2) 



dvi 



, ifl<i,j<N, 

ifl<i<N,j = N + l, 

if i = N + 1,1< j < N, 
ifi = N + l,j = N+l, 

dv2 exp{—vl — {vi — V2)'^} 



■UI+U2 



rU2 /•UI+U2 

/ dvi / dv2 exp{-vl - {vi - ^2)^} 



N 



(2.21) 



(2.22) 



(2.23) 



(2.24) 



l<i<j<N i=l 

Then we can obtain the following result. 

Theorem 2.5 (i) For any fixed z e N;^ and T > 0, as L ^ 00, 'jl^ xi') converges weakly to the law of 
the temporally inhomogeneous diffusion process X(i) = {Xi{t),X2{t), . . . ,Xjv(t)), t € [0,T], with transition 
density gjf{s,x,t,y); 

ffi;(0,0,i,y) = cjvr^VV^(2iv+i)/2^^p|_^|/J^(y)^^(j._^^y)^ (2.25) 
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AfN{T-s,x) 

for < s < i < T, X, y e R+< , where cjv = 1/ H^i r(i). 
(ii) The diffusion process X(t) solves the following equation: 



(2.26) 



Xi{t)=Bi{t)+ I h^{s,X{s))ds, te[0,T], i = l,2,...,7V, (2.27) 
Jo 



where 



d 

b[{t,x) = — lnMv(T-t,x), i = 1,2,. ..,7V. 



Next we consider the case that T = Tj^ goes to infinity as L ^ oo. 

Corollary 2.6 (i) Let Tl be an increase function of L with Tl ^ oo as L oo. For any fixed z G N;^, 
as L ^ oo. fi^ Tl^') converges weakly to the law of the temporally homogeneous diffusion process Y(<:) = 
{Yi{t),Y2{t), . . . , YN{t)), t e [0, oo), with transition density pn{s, x, t, y); 

p^(0,0,^,y) = ?;v^-^(^^+^)/2exp|-^|/i;v(y)^ (2.28) 

Pjv(s,x,t,y) = /jv(i - s,y|x)/ijv(y), (2.29) 

/ijv(x) 

/or < s < i < oo, X, y e R+:^, where c'n = (2/7r)^/V nf^i r(2i). 
(ii) The diffusion process Y{t) solves the following equation: 

Yi{t)=Bi{t)+ [ J—ds+ y If ^-^ ds+ f ^ ^-^ ds\, (2.30) 

Jo Y,{s) i<,.f^,,.^aio F^(s)-y,-(5) Jo Y,{s)+Yj{s) j 

tG [0,oo),z= 1,2,..., TV. 



2.3 Remcirks 

(i) The process X(i) ( resp. X(i) ) represents the system of N Brownian motions ( resp. A'' Brownian 
meanders ) started from the origin conditioned not to collide up to time T. A limit theorem for one- 
dimensional random walk conditioned to stay positive was firstly observed by Spitzer [37] and then studied 
and generalized by many probabilists [1, 12, 3, 6]. For two-dimensional random walk conditioned to stay in 
a cone, a limit theorem was proved by Shimura [36]. Our theorems are multi-dimensional versions of these 
limit theorems. 

(ii) The process Y(i) ( resp. Y{t) ) represents the system of A'' Brownian motions ( resp. A'' three-dimensional 
Bessel processes ) conditioned never to collide. The function /lAr(x) ( resp. h]y{x.) ) is a strictly positive 
harmonic function for the absorbing Brownian motions in the Weyl chamber R;^ (resp. R+;^ ). The process 
Y{t) ( resp. Y{t) ) is the corresponding Doob h-transform [5, 11]. A functional central limit theorem to the 
process Y{t) was also discussed in a recent paper by O'Connell and Yor [31]. 

(iii) The relation between the Brownian meander and the three-dimensional Bessel process was discussed in 
Imhof [13]. From our results Imhof's relation is generalized as follows: 

(Without wall restriction) For any to = < < • • • < = T, ^ e Z+, 

t ^1 

n5?^(*i-i'yi-i'*i'yi) = Civr^('^-i)/4[]pAr(Vi,yi_i,ti,yi)^^--^, (2.31) 
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for any e R< , i = 1,2, . . . ,£, where yo = and 



- CAT N/2r\ 

L. - - 



N 
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L\r(i/2)- 



(With wall restriction) For any to = < ti < ■ ■ ■ < te = T, £ G Z_^_ 

J]5?;(ii_i,yi_i,ti,yi) =?jvT^'/2 JJp;v(ti-i,yi-i,ti,yi):^— -, (2.32) 



for any y^ e R+^, i = 1, 2, . . . , where yo = and 



^- = p; = UJ nro) 



(iv) Consider an ensemble of N x N complex Hermitian matrices {H}. The Gaussian unitary ensemble 
(GUE) is the ensemble with the probability density function 

M^^^(i?,a?)=ciexp|-2^IVF2|, 

where af is variance and ci = 2^^/^(7rcr^)~^^/^. The Gaussian orthogonal ensemble (GOE) is defined as 
the ensemble of N x N real symmetric matrices {A} with the probability density function 

M^°^(A^T2) = ^2exp|-^TrA2| ^ 

where a2 is variance and C2 = 2~''^/2(7rcr|)~''^(^+^)/^. It is known that the distributions of eigenvalues 
X = {xi,X2, ■ ■ ■ ,xn) of these matrix ensembles are given as 

,«^^(x,a?) = |y-^'^^exp{-g|/.;v(x)^ 

and 

GOE/ 2\ Cat _Ar(7V+l)/2 f 1^1 1 . / \ 

5 (x,T2) = -^0-2 exp|-^|/iAr(x)> 
respectively [26]. Theorem 2.2 and Corollary 2.3 give the relation 

5^(0, 0, T, y) = Nlg'^O^iy, T), y e R^, 

and 

PAr(0, 0, t, y) - iV!.g<=^^(y, t), y e R^, t > 0. 

In order to study a Gaussian ensemble of complex Hermitian matrices intermediate between GUE and GOE, 
Pandey and Mehta considered the following probability density functions with a parameter a G [0, 1] 

/.^^(if, a)= JdA 11^"^ {H - A, 2a\2)^GOi5(^^ 2(1 - a^^), 

where = 1/{2(1 + a^)} [27, 33]. They have studied a transition from the GOE to the GUE observed as 
a changes from to 1. Let g^^ {x,a) be the probability density function of eigenvalues in this ensemble of 
Pandey and Mehta. We can show the ec^uality [18] 
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It is shown in [19] that as a consequence of this equality, the Harish-Chandra formula for an integral over 
the unitary group can be obtained. Similar argument concerning the relation between the process X(t) and 

the nonstandard classes of random matrices is given in [20]. 

(v) Spohn [38] constructed nonintersecting Brownian motions on a torus and discussed the infinite volume 
limit to an infinite system of Dyson- type Brownian motions, which was also constructed by Dirichlct form 
technique in Osada [32]. For the present A'' nonintersecting Brownian motion X(i) in a finite time interval 
(0,T], two types of temporally inhomogeneous infinite particle systems are obtained by setting T = T{N) 
and taking N oo. If we set T{N) — 2N and observe the bulk configuration of particles at time t = 
T{N) + s, — oo < s < 0, a spatially homogeneous but temporally inhomogeneous system is derived in the 
infinite particle limit, whose multitime correlation functions have the quaternion determinantal expressions 
with sine-kernel. If we set T{N) = 2N^/^ and the particle configuration at time t = T{N) + s, — oo < s < 
around the position 27V^/^ — is observed, a spatially and temporally inhomogeneous system is derived in 
A'^ ^ 00 , in which multitime correlation functions are given by the quaternion determinants with Airy- kernel 
[30, 17]. It is easier to prove the limit theorems for Dyson's Brownian motion model Y{t) corresponding 
to the above two kinds of infinite particle limits. The former limit provides a homogeneous infinite system, 
which coincides with the system studied by Spohn [38], Osada [32] and Nagao and Forrester [29], and the 
latter does a temporally homogeneous but spatially inhomogeneous infinite system, which is related with the 
process recently studied by Prahofer and Spohn [34] and Johansson [14]. See Nagao [28] for A'' — > oo limit 
of the process X(t). 



3 Proof of Theorems 
3.1 Proof of Theorem 2.2 

Let A^jv(t77, v[u), u,v e Z^, be the total number of the vicious walks, in which the N walkers start from 
Ui,i = 1,2, ... ,N, and arrive at the positions Vi,i = 1,2, ... ,N, at time m. Then the probability that such 
vicious walks with fixed end-points are realized in all possible random walks started from the given initial 
configuration is A''jv(m, v|u)/2'"^, which is denoted by VAr(m, v]u). We also put 

Fjv(m|u) = ^ Fjv(m,v]u). 
vez^ 

Define a subset of the square lattice Z^, 

Cm = {{x,y) e Z'^ : X + y = even, 0<y <m}, 

and £m be the set of all oriented edges which connect the nearest-neighbor pairs {{x,y), {x' ,y')) of vertices 
with y' = y + 1 m Cm- Then each walk of the i-th walker can be represented as a sequence of successive 
edges connecting vertices (uj, 0) and {vi, m) on {Cm, £m), which we call the lattice path running from (wj, 0) 
to {vi,m). If such lattice paths share a common vertex, they are said to intersect. Under the vicious walk 
condition, what we consider is a set of all A?'-tuples of nonintersecting lattice paths. Let 7ro({(Mi, 0)}^]^ — > 
{{vi,m)}^i) be the set of all A^-tuples (tti, . . . of nonintersecting lattice paths, in which tTj runs from 
(wj,0) to {vi,m), i= 1,2, ... ,N. NN{m,v\u) = |7ro({(wi, 0)}^^ {(^i) and the Karlin-McGregor 
formula [15, 16] gives 

iVjv(m,v|u) = ^ det (|7r((uj,0) ^ {vi,m))\\, 

where \A\ denotes the cardinality of a set A and 7r((ttj,0) {vi,m)) the set of lattice paths from (ttj,0) 
to (vi,m). (Such a determinantal formula is also known as the Lindstrom-Gessel-Viennot formula in the 
enumerative combinatorics, see [23, 10, 39].) Since \Tr{{uj,0) {vi,m))\ = {(^m+u"-v)/2) ' have the 
binomial determinant 

yiv(m,v]u) =2-™^ det ((^ ^ V (3-1) 
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For L > we introduce the following functions: 

Lx 



, xeR, and (j)L{x) = {(l)L{xi),(j)L{x2), ... ,(l)LixN)), xeR 



N 



2 

where [a] denotes the largest integer not greater than a. We show the following lemmas. 
Lemma 3.1 (i) Fort>0,yi€ and y G 

(§)'^^iv(0L^(t),0L(y)|x) =4t-^^/2/z^(^)exp|-^|/.^(y) (^l + O (^M^) , (3.2) 

OS L ^ 00, where c'j^ ^ {'^'^)~^^'^/Uf=i^{j)- 
(ii) Fort>0 and x e 

VN{(I^Lm^) = ^flN ( (l + 0(y)), (3.3) 



as 00, where cn = tt^^' nf=i{r(j)/r072)}. 

Proof. It is enough to consider the case that x = 2u, (pL{y) = 2v, u, v e Z^ and ^^,2 (t) = 2£, i G Then 
Nr,{^L<t),My)\^) = Afiv(2^,2v|2u) = det (( 

and 

2^ \ _ {2ey. 



£ + Uj -ViJ {£ + Uj - Viy.{£ - Uj + Uj)! 

(2£V 

with 

where (a)o = 1, (a)fe = 0(0 + 1) • • • (a + — 1), > 1. Then 

JX- C2/'V 

The leading term of deti<jj<Ar(>lij(£, v, u)) in L ^ 00 is 

'i + v, 



Di{v,u) = dot 



l<i,j<N \ \£ — V, 



(_l)W(A^-i)/2 det f f 

l<i,j<JV \ \ 



£ + V^ 



l<i,j<N W£ — Vi^ 

Let ^(u) = (^i(u), . . . , ^Ar(u)) be a partition specified by the starting point 2u defined by 

^,(u) = UN-j+i - (iV - j), J = 1, 2, . . . , iV. (3.5) 
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Noting that the Vandermonde determinant deti<j,j<jv(2;i^ "') = ni<i<i<]v(^i ~ ^j)> have 

l<i<j<N ^ * J / ^ ^ -"^ / 

n 2e{vj - Vj) f e + vi e + vrA 

i<i<j<N - ^^)(^ - ''^"^ - ^1 ' • • • ' ^ - t;;v j ' 
where sa(^1) ■ ■ ■ , zn) is the Schur function associated to a partition A = (Ai, A2, . . . , \n) defined by 

deti<ij<]v {z^'"^'^~'j 

sx{zi,...,zn) = / ■ (3-6) 

deti<i,j<jv Uj •'j 

(See Macdonald [25].) It is a symmetric polynomial of degree -zi, • ■ • , and it is known that 

(see p.44 in [25]) 

aA(l,l,...,l)= n (3.7) 

l<i<j<N ■' 



Then the leading term of £>i(v, u) in L — » 00 is 



l<i<j<JV 



= r^(^-l)/^2^(^-l)/2/l;v(v)/l;v(u) n 



7 — % 

l<i<j<N 

hN{j)hN{2u)f[^y (3.8) 



By Stirling's formula we see that 



n,.^— "'n(>-S)""'(^)"0-G))- - 

From (3.4), (3.8) and (3.9) 

FAr((/>L2(t),0i(y)|x) = 2-2^^iVjv(0L2(t),,^i(y)|x) 

-w(^)"^«(^)*~(-)»=p{-^}(i.o(M)) 

= c'„(|)''r-/.,„g)e.,{_M!},.„(„(,,,(M); 



Then we obtain (3.2). 

By (3.2) and simple calculation we have 



ViMt)\^) = J^^ e-\y\'/'%N{y)\ (1 + O , (3.10) 
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as L — > oo. The last integral is the special case (7 = 1/2 and a = l/2t) of 

/ rfue-«H^ TT 1^ - H,|27 = (2^)iv/2(2^)-Jv(T,(iv-i)+i)/2 tt r(l + 11) ^3_^^^ 

found in Mehta (eq.(17.6.7) on page 354 in [26]), whose proof was given in [24]. Then we have (3.3) by 
elementary calculation. This completes the proof. | 

Lemma 3.2 Let f > and x, y e R^. Then 

VN{Mt),H{y)\H{^)) = fN{t,y\^)(l + o{^^^^\], (3.12) 



as L ^ 00. 
Proof. From (3.1) 



2] VN{Mt)AL{y)\<l>L{^)) 
2-^4>.^it)(L\'' ff Mt) 



det ( 2-*i-=W-iL 



1<IT<N y \{(pL^it) + (pL{Xj) - <l)L{yi))/2 

Application of Stirling's formula yields the lemma. | 

By Donskcr's theorem (see, for instance, Billingsley [2]) wc sec that J\f]\f{t,x.) is the probability that 
Brownian motions do not collide until time t. We have the following asymptotic behaviours of the function 
jVjv(t,x) as |x|/\/t 0. 

Lemma 3.3 Lett>0 and x e R< . Then 

^„(,,, = _L,„(^)(i,e.(W)). W^o. (3.13) 

wkere c„ = n?=i{rO)/r072)}. 

Proof. First note that 

/;v(i,y|x) = (27rt)-^/2exp|-l^(a.,2 +t/,2)| ^^det^^ (e^^^'*/*) . 

We rewrite the determinant as 

dot fe^^-/*) = det:<,,-<. iie^^/y^-^^^) ^ (r,..fuN-A 
i<t,j<N V / deti<,j<./v {{e"='/^}^-3) i<i,j<N V ) 

= s^(y) (e"^'",e"^/',...,e"-^'/') ]J (e^^V* - e"'/*), 

l<i<j<jV 

where ^i(y) = yN-i+i - {N - i), i = 1,2, . . . ,N. Using it 

fN{t, yjx) = {2nt)-^^h^^y) (e-^/*, e^^/*, . . . , e^"/*) 

|2 I |,J2 



xexp<|- ' ' |> [] (e-.V* - e^^/*). (3.14) 

i<j<j<Ar 
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Since 

N 



lim S5(y)(e-/*, . . . , e--/*) = s^^^^il, 1, . . . , 1) = h^iy) f[ -i-, 



and 



l<i<j<Ar 
the function is asymptotically 



n (e^^/*-e-/*)=/..(f)(l + o(M)), M^o, 



4..M.)exp{-f}(...(^)), 



By (3.11) we have (3.13). | 

Lemma 3.4 For z e andT>0 {/xf, i > 1} is tight. 

Proof. By the Kolmogorov's tightness criterion it is enough to prove that for any s > 



limsupA<i,T ( ^ max \w{u) - w{v)\ > e ] = 0. (3.15) 

(See, for example, Billingsley [2].) Since /xf^ rp is the probability measure of a linearly interpolated random 
process, (3.15) is derived from the following estimates: as 5 ^ 0, 

limsup/x^ T ( max \w{u) — w{0)\ > e/2 ) = o{S), (3.16) 

I,_>oo ' \0<u<S J 

limsupAi^r ( max |«;(t + w)-M;(t)| >e/2 ) =o(5), t€[d/2,T-S]. (3.17) 

L — oo ' \0<ti<5 J 

Under the nonintersecting condition, for any i ~ 1,2, . . . , N 

\wi{u) - Wi{0)\ < \wn{0) - wi{0)\ + {wi{u) - wi{0))- + {wn{u) - wn{0))+, (3.18) 
where a+ = max{a,0} and a_ = max{— a, 0}. Then the set 

max \w{u) — w{0)\ > s/2 

0<u<S 

is included in the set 

< max (wi(u) — wi(0))- > — - \- U \ max (wn(u) — WAr(0))+ > — - > . 

Noting that {wi{u) — 'u;i(0))_ and {wn{u) — WAr(0))+ are nonnegative submartingales, we can apply Doob's 
theorem (see, for instance, Revus and Yor [35]) to obtain 

< (^^y ElT{\wi{S)-wim'' + \wN{S)-WN{OW) (3-19) 
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for any p > 1 and L > A{zn — z\)le, where E\ j, represents the expectation with respect to the probabiUty 
measure /i^ p. 

From Lemmas 3.1, 3.2 and 3.3 



Umsupi;^ rfkiW - ^"i(0)|^ + \wn{S) - WN{^)f\ 



< C2C^5-^^/^ J^^ dy exp {-^1 hNiyfiy^, + vl) 

< Cs5P'^ j^^ dx exp {-^} hN{^f{x{ + 
= 0{SP/^). 

Taking p > 2, we obtain (3.16). 

Fix t e [S/2, T-6]. By the Markov property 

l^l,T ( max^ \w{t + w{t)\ > 

SiL\t + u))- S{LH) 



= Qhr \ max 



- P^(Ai2T) \^ ^ «' \0<u<5 



> e/2 
S{L'^u) - 5(0) 



> £/2 



(3.20) 



where = {Sj < Sj < ■ ■ ■ < Sj' , < j < m} and represents the expectation with respect to the 
probability measure P^. By Doob's inequahty for any x G and p > 1 



P~ I max 

. 0<u<5 



S{L^u) - S{0) 



>£/2 ) < ( '^Y 



S{L^d) - 5(0) 



(3.21) 



Prom Lemma 3.1 (ii) 



Umsup/i^ rp I max \w{t + u) — w{t)\ > e/2 ) 

L^oo ' \0<u<5 J 

L^oo P"(Ai2r) ^ ^ 
Taking p > N{N - l)/2 + 2, we obtain (3.17). This completes the proof. | 



(3.22) 



Lemma 3.5 Let z G Z'^, = to < h < ■ ■ ■ < tk-i <tk=T and9 = (Oi, ...,0k)e K'^". Then 



lim Ef^ y exp < V^y^ 6j ■ w{tj 



-i 



fe 



dyidy2 ■■■dykT[ gN{tj-i,yj-i,tj,yj) exp <^ V Oj ■ yj 



(3.23) 



where yo = 0. 
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Proof. By Lemmas 3.2 and 3.1 



= lim 



1 



^L^{t2-ti),S{(j)L2{t2 -ti)) =X2,- 



= lim 



1 



^L^{tk-tk-i)^ S{(j)L-^{tk - tk-i) = Xfc,exp ^ • 



r — pz^A ^ / dyidy2---dyk , 



Affc 



VN{Mti)AL{yi)\^) 



xVN{<l>L^{t2 - ti),<l)L{y2ML{yi)) ■ ■ ■ 
xVNi(f)L^{tk - ife-l),0L(yfe)|0L(yfe-l))exp I y/^^0j ■ yj 



i=i 



= / dyidy2 ■■■dykY] Yj-i, Yj) exp < V^Y] 0, ■ y. 

This completes the proof. | 

Proof of Theorem 2.2. Prom Lemmas 3.4 and 3.5 we see that fJ-L.ri') converges weakly to the law of the 

time inhomogencous diffusion process X(i) with transition density g^(s,x, t.y). Noting that J^f^iT — t,x) 
is a solution of the heat equation, we see that gjf{s,x,t,y) satisfies the following backward equation: 



d 1 

—gjf{s,x,t,y) = -A^g%{s,yi,t,y) + b'^{t,x) ■ V^g%{s,x,t,y). 
Then the process X(t) solves (2.12). This completes the proof of Theorem 2.2. 



(3.24) 



Proof of Corollary 2.3. By simple observation we see that the estimates concerning the probability ii\ j, 
in the lemmas are uniform with respects to T. Then (i) is obtained from the properties 



lim gj^{0, 0, t, y) = pjv(0, 0, t, y), 

1 — xyo 



lim g^{s,x,t,y) ^ pNis,x,t,y), 

T — *oc 



which are derived from Lemma 3.3. Noting that /ijv(x) is a harmonic function, we see that pjv(s, x, y) 
satisfies the following backward equation: 

d 1 

—pN{s,x,t,y) = -AxPjv(s,x,i,y) + Vxln/ijv(x) • VxPjv(s,x,t,y). (3.25) 
Then the process Y{t) solves (2.15). This completes the proof of Corollary 2.3. | 

3.2 Proof of Theorem 2.5 

Let iVjv(m, v|u), u, v e N^, be the total number of the vicious walks with wall restriction, in which the A'' 
walkers start form the positions u,, i = 1, 2, . . . , A^, and arrive at the positions Vi,i= 1,2, ... ,N, at time m. 
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Then the probabiUty that such vicious walks with fixed end-points are reahzed in all possible random walks 
started from the given initial configuration is Nj^{m,w\vL)/2'^^ , which is denoted by VatI^ti, v|u). We also 
put 

VAr(m|u) = ^ yAr(m,v|u). 
By the Karlin-McGregor (Lindstrom-Gessel-Viennot) formula, we have [22] 



Vjv(m,v|u) = 2-"''^ det , x ,„ - L x J • (3-26) 

^ ' ^ i<i,j<N\\{m + Uj -Vi)/2j \{m + Uj+Vi)/2 + lJ J ^ ' 

Let 

for a partition A = (Ai, . . . , Xn)- Remark that spx{zi, . . . , Zn) is the character of the irreducible represen- 
tation corresponding to a partition A of the symplectic Lie algebra (see, for example, Lectures 6 and 24 in 
Fulton and Harris [9]). By using the function spx instead of the Schur function, we can show the following 
lemma by a similar way to the proof of Lemma 3.1. 



Lemma 3.6 (i) For t > 0, x e and y e R+< 

[^yVM{<l>L<t),My)\^) = c'ivi-^(^^+^)/^/i^ (j) exp|-^| ft^(y) (^1 + (^M^^ , (3.28) 



as 00, where c'n = (2/7r)^/V nV=i r(2j)- 
(u) For i > and x e N 



N 



M0.^W|x) = i-/i.(^)(l + o(l)), (3.29) 
as 00, where cn = {7r/2)^/^U^^,{T{2j)/r{j)}. 

Proof. Again we consider the case that x = 2u, <l)L{y) = 2v, u, v G and (pL^it) = 2£, £ G Z^. By the 
equation (3.26) 



AT 

n {£ - Viy.{e + vi)! i<Tj 

with 



Then the leading term of deti<,j<jv(^ij(^, v, u)) as L — > 00 is 

^ ^ i<i,j<N \ \£-ViJ \e + V: 



^ ' l<i,3<N \\£ - Vi J \£ + tHj J 



= det 



^<i,0<N \\l + V^ \i-V^J 1 '^^^'''> - Vl ' t - Vn 
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where^(u) = (^(u), . . . ,CAr(u)) with ^j{u) = un-j+i - {N - j + = 1,2, . . . , N. Note that 

7=1 ^ ■'^ l<i<7"<Af \ » J / 



N 

l<i,j<N 

Then by simple calculation we have 



It is known that 



det , , , , , „ 

l<i,j<N \\£ + Vi J \i-Vi 



ap,(i,...,i)= n 4-An- 



with = Aj- + iV - j + 1, mj = A'' - j + 1 [9]. Then 

l<i<i<N ■' 1=1 ,=1 w ' 



.,2 „2 AT AT 



l<i<j<N •' j 

Then the leading term of Di{v,u) in L — > oo is 



^2(V,U) = /I;v(u)n— Y II 



^ 2 



=nfM'"(|)^-(2v)- 



Prom (3.9), (3.30) and (3.34) we have 

^ 1 / ^ \ 



= (I) " (i) S„w e.p {-f } (: . o (M)) 



Then we obtain (3.28). 

By (3.28) and simple calculation wc have 



P 4--Ar(2Ar+l)/2 / y \ /" / /l 

n^vWIx) = '-^^y^'^- (I) X„ .-y e-" '»fc(y)| (l + O (i 

= r(lTT)"*»(z^)t''--'-'"'^|S"«l('-Kl))' 



as L — > oo. The last integral is the special case (7 = 1/2 and a = 1) of 

r ^ 

/ due-HV2 -Q |„2_^2|2.-Q 
•'^'^ l<i<j<N j = l 



N 

L, |2o-l 



N 

r(i + 7) 



r(l+j7)r(a + 7(j-l)) 



J=l 
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((17.6.6) on page 354 in [26]). Then we obtain (3.29) by elementary calculation. | 

Following the same calculation as was done in the proof of Lemma 3.2, we have the following lemma. 

Lemma 3.7 Let i > and x, y e R+< • Then 

Vn (0L.(t),0L(y)|0L(x)) = h{t,y\^) {l + O (^^^)) , (3.36) 
as L —> DO. 

We rewrite /Ar(i,y|x) as 

7^(t,y|x) = (27r0-^/^sp,-(,) (e--/*, . . . , e^^/*) exp | - + | 

X J|(e^^/*-e-^^/*) I (e^i/*-e^'/*)(e(^'+^^)/*-l) i J [|e=^^/* 

j=l l_l<i<j<JV J [j=l 

Then wc can obtain the following lemma by a similar way to prove Lemma 3.3 by virtue of the equations 
(3.32) and (3.35). 

Lemma 3.8 Let f > and x e R+< . Then 

M^{t,^) = lSh^(2L\(i + o(^\\, M^o, (3.37) 



-N+l 



cn \ViJ \ \ViJ J ' Vi 

where = (7r/2)^/2 nf=i{r(2j)/r(i)}. 

From the above lemmas and the same argument as in the previous subsection we can obtain Theorem 2.5 
and Corollary 2.6. 

4 Proof of Lemmas 2.1 and 2.4 

We use the following identity, which is shown in de Bruijn [4]. (See also Appendix in [18].) Lemmas 2.1 and 
2.4 are easy consequences of this result as shown below. 

Lemma 4.1 Let z be a square integrable piecewise continuous function on R^. Then 

[ dy det (^(a;,,y,))=Pfi<.<.<^(F,,(x)), (4.1) 

where 

jCr _ ( N, if N is even 

\ N + l, ifN is odd, 

/OO 
z{xi,y)dy, 
-OO 

I rfT"- lt■^\ ■rf'r- \ 

dy\dyi, 
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I,{x,,Xj), if l<i<j<N, 

-h{xi,Xj), if 1 < j < i < N, 

Izixi), if l<^<iV,J=A^+l, 

-h{x,), if i = N + l,l<j <N, 

0, if l<i=j<N+l. 



Proof of Lemma 2.1. From the above lemma and integration by substitution, it is enough to show 

hAxi) = h hAxi,Xi) = '^(^^^y (4-2) 
for zi{x^y) = e~^^~^^^/\/7r. The first equation in (4.2) is trivial. Let zo{x) = e~^^ /y/ir. Then 

/OO nXj—Xi 
dyi dy2 zo{yi)zo{yi+y2) 

-OO Jxi—Xj 
-t fXj—Xi pOO 

= - dy2 e-«i/2 / ^-2iy,+y,/2r 

^ Jxj—X^ J — OO 



1 pXj—Xi / ^ 

/ dy2 e-^^'/2 = * ( - 
V27r Jxi-xi \ 



V2 



This completes the proof. | 

Proof of Lemma 2.4- It is enough to show 



Iz2{xi) = {xi) , h^{xi,Xj) = -^5/ {xi,Xj) , (4.3) 

for Z2{x,y) = K^e-^^-?/)" - e-^^+J')') /y^j l{x > 0,y > 0}. The first equation in (4.3) is obtained easily. 
To show the second equation we put 

G((ai,a2],(6i,62]) = / dyi I dy2 zo{yi)zo{yi - y2), 

for ai, a2, bi, 62 € R U {—00, 00}. Then 

Iz2{xi,Xj) = G{{-Xi, 00), {-00, Xj - Xi]) - G{{-Xi, 00), {-00, -Xj - Xi]) 

— G{{xi, 00), {-00, Xj +Xi]) + G{{xi, 00), {-00, -Xj +Xi]) 

— G{{~Xj, 00), (—00, Xi — Xj\) + G{{—Xj, 00), (—00, —Xj — Xi]) 
+ G((a:j, 00), (-00, a;^ + Xi]) - G{{xj,oo), (-00, -Xi + Xj]) 

— G(( Xi , Xi^ , (^Xi Xj , Xj '^i] ) G(^(^Xi , Xj^ , i^Xj Xi , Xi -\- ] ) 

G( ( Xj , Xj] , ( Xi Xj , ] ) 

This completes the proof. | 
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